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ABSTRACT 

We present results for a suite of fourteen three-dimensional, high resolution hydrodynami- 
cal simulations of delayed-detonation models of Type la supernova (SN la) explosions. This 
model suite comprises the first set of three-dimensional SN la simulations with detailed iso- 
topic yield information. As such, it may serve as a database for Chandrasekhar-mass delayed- 
detonation model nucleosynthetic yields and for deriving synthetic observables such as spec- 
tra and light curves. We employ a physically motivated, stochastic model based on turbulent 
velocity fluctuations and fuel density to calculate in situ the deflagration to detonation transi- 
tion (DDT) probabilities. To obtain different strengths of the deflagration phase and thereby 
different degrees of pre-expansion, we have chosen a sequence of initial models with 1,3, 
5,10, 20, 40, 100, 150, 200, 300, and 1600 (two different realizations) ignition kernels in a 
hydrostatic white dwarf with central density of 2.9 x 10^ g cm~^, plus in addition one high 
central density (5.5 x 10^ g cm~^) and one low central density (1.0 x 10^ g cm~^) rendition 
of the 100 ignition kernel configuration. For each simulation we determined detailed nucle- 
osynthetic yields by post-processing 10^ tracer particles with a 384 nuclide reaction network. 
All delayed detonation models result in explosions unbinding the white dwarf, producing a 
range of ^^Ni masses from 0.32 to 1.11 Mq. As a general trend, the models predict that the 
stable neutron-rich iron group isotopes are not found at the lowest velocities, but rather at 
intermediate velocities (~3,000— 10,000 km s~^) in a shell surrounding a ^^Ni-rich core. 
The models further predict relatively low velocity oxygen and carbon, with typical minimum 
velocities around 4,000 and 10,000 km s~^ respectively. 

Key words: nuclear reactions, nucleosynthesis, abundances — supernovae: general — white 
dwarfs 



1 INTRODUCTION 

Type la supernovae (SNe la) play essential roles in the basic 
frameworks of many branches of astrophysics: In star forma- 
tion and galaxy dynam ics by heating cold interstellar gas (e.g. 
IScannapieco et al.ll2008h, in high energy astrophysi c s as sources of 
Galac tic positrons (e.g. IChan & LingenfelteJll993l : IPrantzos et alj 
[mil), in galactic chemical evolution by enriching the interstel- 
lar gas with a, F e-peak, and possibly p-process eleme nts (e.g. 
Timmes et alj 1 1995: Kobayashi & Nomoto 2009; Travagl io et alj 



20111) and last but not least in cosmology as distance indicators 



re.g. lRiess et alJI 19981 : 1 Schmidt et al.lll998l : IPerlmutter et al.lll999h . 
In spite of their ubiquitous presence in astrophysics, no progen- 
itor systems have been observed and unambiguous identification 
of their nature remains elusive. Essentially by means of exclusion, 
white dwarf (WD) stars in interacting binary systenis are the only 
viable proposed progenitor systems (cf. lBloom et al.|[2012b . 

In the last few years we have witnessed revived in- 
terest of the double detonation He-accretion channel (e.g . 
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ing support from theorists for the double degenerate channel 
dOilfanov & Bogdanl201Ql : IPakmor et al.l201Ql : lpSmor et al.l201ll : 
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IPakmor et alJ l2012b . culminating in mounting evidence that SN 
201 Ife was possib l y due to a merger of two white dwarfs (e.g. 
iBloom etal.ll2012l : Ichomiuk et aP l2012l : iRopke et all l2012h . In 
spite of these recent developments, Chandrasekhar-mass models, 
which had been the favoured explosion scenario by many in the last 
two decades, still retain a list of strong arguments in their favor. 
For example, blue shifted Na-absorption features along the lines 
of sight towards SNe la a re interpreted as a clear signature of th e 
single degenerate channel (jPatat et al.ll2007l : [Sternberg etaDl201lh . 
Further support to the single degenerate scenario is given by the 
fact that two recurrent nova systems are known where the accreting 
WD is near the Chandrasekhar-limit (RS Oph and U Sco), which 
tells us that potential progenitor systems do exist in nature. In fact, 
the supernova P TF 1 Ikx is best de scribed by a symbiotic nova pro- 
genitor system (iDildav et al.ll2012h . 

The issue is further complicated by the fact that the evolu- 
tion of zero-age binary systems towards potential SN la progenitor 
systems is still not well understood and remains a very active area 
of research today (e.g. iRuiter et al ] l2009l . I2OIII : iMennekens et"aD 



fined detonation (GCD) Jjordan et"aD 



2008lJ2012h. and pulsational 



l20ld : IWang^rani20K)h . Binary population simulations, in which 
a large number (>10^) of binaries can be rapidly evolved from 
the zero-age main sequence for a Hubble time, are the only way 
in which one can obtain reliable estimates of relative birthrates 
for different progenitor scenarios. However, the physics of binary 
star evolution is very co mplex and some evolution ary phases are 
poorly understood (e.g.. Ivan der Sluvs et aLllioiol) . For example, 
the evolution of single degenerate Chandrasekhar mass progenitors 
strongly relies on the assumptions made about WD mass accretion 
rates and retentio n efficiencies, for which there are differing camps 
of th ought (cf. IPrialnik & Kovetzl 1 19951 : [ Han & Podsiadlowskil 
I2OO4I) . Despite the different assumptions made in various binary 
population codes, many of the codes do predict that single degen- 
erate (Chandrasekhar mass) scenario progenitors are still promis- 
ing candidates for at l east some fraction of Type la supernovae (see 
iNelemans et"aDl20ll table 1). 

Near Chandrasekhar-mass explosion models in the single de- 
generate channel have long been considered as favorites to ex- 
plain SNe la. The realization that a detonation burning through a 
near Chandrasekhar-mass white dwarf in hydrostatic equilibrium 
jArnett 1969) produces mainly material that has been processed to 
nuclear statistical equilibrium (NSE) and not enough intermediate 
mass elements (IMEs), such as silicon or su lfur, has lead to t he in- 
troduction of "delayed detonation" models (lKhokhlovlll98^ . The 
key features of delayed detonation models are the following: 



• First nuclear burning is ignited in a deflagration flame produc- 
ing mainly iron group elements (IGEs) in the initial burning phase 
at high density. 

• The energy released in this sub-sonically propagating mode of 
nuclear burning leads to an expansion of the star, moving unburned 
nuclear fuel to lower density. 

• After some time delay a supersonically moving mode of nu- 
clear burning - a detonation - emerges. 

• The supersonically moving detonation front quickly overruns 
much of the remaining fuel, a significant fraction of which only 
burns to IMEs owing to the reduced burning densities resulting 
from the pre-expansion. 

Several three-dimensional models of delayed detonation 
Chandrasekhar-mass explosions have been publish ed in the vari- 
ants o f deflagratio n to detonation t r ansition (DDT) jGamezo et al 
20051: Ik opke & Niemeve j l2007l: Isravo & Garcia-Senj 12001 ' 



Seitenzahl et alJ 1201 ll : iRopke et all l2012h ~ gravitationally con 



reverse detonatio n (PRD) models (iBravo & Garcia- Si^ l2009l : 
iBravo et aDl2009l) . Generally, information on the chemical com- 
position of the ejecta for these models is either not given or limited 
to ^^Ni and a coarse description of elemental yields or major nu- 
cleosynthesis groups, s uch as unburned car bon or oxygen, IMEs, or 
IGEs. The exception are lBravo et al.1 (l2009l) . who show a table of 24 
isotopes for two PRD explosion models. While full isotopic infor- 
mation of the ejecta for three-di mensional pure deflagration explo- 
sions e xists in the literature fe.g. lTravaglio et al ]|2004l : lRopke et alJ 
l2006ah such detailed information for delayed detonation explo- 
sion models is currently only available fo r a few two-dimensional 
(iMeakin et al.ll2009l : iMaeda et a l. 2010b) or one-dimensional ex - 
plosion models (e.g. Ilwamoto et al. 19991 : iBrachwitz et"aDl200Qh . 
Here, for the first time, we present detailed nucleosynthetic yields 
for a suite of high resolution, three dimensional hydrodynamical 
explosion simulations, producing a large range of ^^Ni masses be- 
tween -0.32 and 1.11 M0. 

In Section [2] we introduce our initial stellar models and the 
ignition setups, briefly describe our thermonuclear supernova hy- 
drodynamics code Leafs, elaborate on how we model the defla- 
gration to detonation transition, and discuss the morphologies of 
the explosion models. In Section [3] we summarize how we obtain 
the full compositional information by post-processing of tracer par- 
ticles and we present the yields for all models. In Section |4l we 
discuss the relevance of our contribution and conclude with an out- 
look. 



2 HYDRODYNAMIC SIMULATIONS 
2.1 Initial models 

For initial stellar models we have chosen isothermal non-rotating 
WDs in hydrostatic equilibrium. The models have a central density 
of pc = 2.9 X 10^ g cm-3 (M = 1.400 M©, R = 1.96 x 10^ cm), 
but we also include a low central density Pc = 1.0 x 10^ g cm~^ 
(M = 1.361 M0, R=1.96x 10^ cm) and a high central density 
Pc = 5.5 X 10^ g cm-^ (M = 1.416 Mq,R = 1.96x 10^ cm) version 
of the NlOO model IRopke et aLt2012) . These low and high central 
density models are labeled NIOOL and NIOOH, respectively. 

For all models we set up a cold (7 = 5x10^ K), hydrostatic 
WD for a homogeneous composition of ^^C and ^^O in equal parts 
by mass. To account for an assumed solar metallicity of the zero- 
age main- sequence progenitor, we set the electron fraction to Yq = 
0.49886, which corresponds to a ^^Ne content of 2.5 per cent by 
mass (see Sec[3])- 

The ignition of the deflagration is critical for the outcome of 
the explosion. This stage, however, is difficult to model. A cen- 
tury of convective carbon burning precedes the thermonuclear run- 
away. The flow field in this "simmering phase" is expected to 
be highly turbulent and thus hard to address numerically (bu t see 
e.g.,|Hoflich & Stein 2002; Kuhlen et al. 2006; Zingale et al. 200i; 
iNonaka et al.lEoi2h and not finally settled. Not explicitly modeling 
the ignition phase, we therefore take a pragmatic approach and treat 
the ignition geometry as a free parameter. This has two aspects: (i) 
the ignition shape has to seed unstable modes in a realistic way 
and (ii) it should capture the correct position of the ignition re- 
gion relative to the WD's center. Regarding (i), we chose to ignite 
the deflagration in a number of spherical sparks near the center of 
the WD. Although recent simulations suggest that a high number 
of nearly simultaneous ignitions appears unlikely (iNonaka et al.l 
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Figure 1. Ignition geometries of the deflagration for all models. We show the arranged ignition kernels (red spheres) and a transparent blue contour where the 
distance to the center is 2.5 x 10^ cm. For the models with pc = 2.9 x 10^ g cm~-^ this radius corresponds to a density of pfuei = 2.2 x 10^ g cm~^. Models 
N300C and N1600C have a very compact and dense arrangement of the ignition kernels resulting in a setup of high spherical symmetry. 
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(e) N3; t = 1.15 s (f) NlOO; t = 1.00 s 



Figure 2. Snapshots of the hydrodynamic evolution of the N3 model (left column) and the NlOO model (right column). The top row shows the rising plumes 
of the deflagration level set (white) during the Rayleigh-Taylor unstable stage of the deflagration phase embedded in a volume rendering of the density (in 
g cm~^). The middle row shows the density and deflagration level set structure at the time the first DDT occurs. The bottom row shows the subsequent 
spreading of the detonation level set (blue) from the DDT initiation sites. 
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l2012h . it has been argued before that such a multi-spot ignition is a 
probable outcome of the convective carbo n simmering phase lead- 
ing up to the therm onuclear runa way (e.g. Garcia- Senz & WooslevI 
I1995I : IWooslevet al . 2004; lapi chino et al.ll2006h . Such a config- 
uration has numerical advantages because it gives rise to a broad 
spectrum of perturbations that seed the Rayleigh-Taylor instabil- 
ity, while smoother ignition geometries leave this to numerical dis - 
cretization errors which depend on resolution (iRopke et alJl2007bh . 
Regarding (ii), in cases with many ignition sparks we set off burn- 
ing on different sides of the center. This corresponds to a nearly 
central ig nition of the deflagra tion. In contrast, the rece nt simu- 
lations of IZingaleetaD (l2009h and iNonaka et all (l2012h suggest 
off-center ignitions due to a dominantly bipolar flow through the 
center. However, even slight rotation would disrupt the bipolar 
flow and consequen tly asymmetric ignition becomes less likely 
(iKuhlen et al.ll2006h . Even in the non-rotating case, for realistic 
Reynolds numbers (~10^^ times larger than what is achieved in 
numerical simulations), the dipole pattern would likely give way to 
a highly turbulent and chaotic flow, which again would suggest a 
more symmetric ignition. 

For this work we investigate a set of fourteen different ex- 
plosion models, corresponding to different ignition setups that are 
summarized in Table [T] These setups are generated from a Monte- 
Carlo based algorithm. The main parameter of the ignition scenario 
is the number of spherical ignition kernels A/^ that are aligned in the 
central area of the WD following a Gaussian distribution in radius. 
To avoid extreme outliers, we require that the ignition kernels have 
to be contained within a sphere of radius R = 2.5g, where a is the 
variance of the Gaussian distribution. We mention that for ignition 
models with low A^^ the placement of the kernels cannot be consid- 
ered as a real Gaussian distribution anymore. In particular model 
Nl with = 1 constitutes a simple single off-center configuration 
where R = G holds. The ignition models N300C and N1600C (rep- 
resenting "compact" ignition configurations) are b ased on a config- 
uration that has been used in lRopke et alj (12007 ah . By the genera- 
tion of these models R = 2.5a holds, but after all ignition kernels 
are placed, the ones with a larger distance than a are removed and 
replaced until for the whole configuration R = a holds (as in the 
case for A^k = 1). As a result we obtain a very dense and compact 
arrangement of the ignition kernels, where the whole configuration 
exhibits a highly spherical symmetry. We use a uniform radius for 
all ignition kernels set to = 10^ cm for all models except for the 
most centrally concentrated and compact ignition model N300C, 
where = 5 x 10^ cm was used. The length defines a mini- 
mum separation distance between the centers of the ignition ker- 
nels. Note that for large A^k, > dy^, hence the ignition kernels 
may partially overlap, in particular in the vicinity of the center of 
the white dwarf. For a visualization of ignition setups see Fig.[T] We 
have chosen the ignition setups to obtain a suite of models covering 
a large range of strengths of the deflagration phase and associated 
pre-expansion of the WD prior to the onset of the detonation. This 
way, we obtain events that cover t he same range of ^^Ni masses 
that is derived for normal SNe la in lStritzinger et al.l (l2006h . 

2.2 Computational method 

All hydrodynamic explosion simulations presented here were per- 
formed with our thermonuclear supernova code LEAFS, a three- 
dimensional finite-volume discretization of the reactive Euler equa- 
tions. The hydro dynamics solver is based on the PROMETHEUS 
implementation jFrvxell et al.l 1989|) of t he "p iecewise parabolic 
method" (PPM) bv lColella & Woodward! (Il984l) . 



Table 1. Parameters of the ignition setup of the deflagration. 
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Deflagration and detonation fronts are modeled as separate 
discontinuities between carbon-oxygen fuel and nuclear ash and 
their propagation is tracked with a level- set scheme ('Rein ecke et alJ 
1999; Osher & Sethian 1988; Smiljanovski et al. 1997). All mate- 
rial crossed by these fronts is converted to nuclear ash with a com- 
position and energy release depending on fuel density. Although 
modeled with the same method, the propagation velocity, the ash 
composition and the nuclear energy release is different for defla- 
grations and detonations at a given fuel densit y. For detonations , 
the corresponding data is taken from the tables o f lFinketal.l(l201Qh . 
and another table for deflagrations was determined in the same way 
as described there. In our hydrodynamic simulations, we model the 
composition of the material with only five species: carbon, oxy- 
gen, a representative proxy for an intermediate-mass isotope, and 
a mixture of nickel and alpha particles representing nuclear statis- 
tical equilibrium configurations. The latter are adjusted according 
to the prevailing thermodynamic conditions. We keep track of neu- 
tronization (and its effect on the equation of state) by following the 

evolution of Yq as a separate and independent passive scalar. 

The speed of the detonations is modeled as in iFink et all 

tOld) : at high densitie s (p > lO'^ g cm__, pathological case), 
speeds are taken from iGamezo et al.l (Il999h : at low densities, 
Chapman-Jouguet like speeds were calculated for the incomplete 
burning yields in our detonation tables. After a very short phase 
of laminar burning following ignition, the propagation of deflagra- 
tions is dominated by buoyancy and shear-induced instabilities and 
interactions with a complex turbulent flow field. The unresolved 
acceleration of the fl ame due to turjjulenc e is accounted for by a 
sub-grid scale model JSchmidt et al.ll2006al lbl). 

To follow the explosion ejecta to 100 s, where homologous 
expansion becomes a good approximation, we employ a moving 
mesh technique (Ropke 2005; Ropke et al. 2006b) with an outer 
coarse inhomogeneous grid tracking the WD's expansion and an in- 
ner finer homogeneous one encompassing the deflagration region. 
We include self-gravity with a monopole gravity solver. The total 
grid resolution is 512x512x512 cells for all simulations presented 
here. 



2.3 Criterion for deflagration to detonation transition (DDT) 

Although details of the transition process from a subsonic 
deflagration to a supersonic detonation are not well under- 
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stood, it is generally believed that strong turbulent veloc- 
ity fluctuations at the deflagration flame a r e required for 



a DDT to occur ( Niemever & Wooslev* Jl997l : iLisewski et alJ 
I2OOQI : IWooslevI laOOvi rfwooslev et al. 2Q09h. Under this condi- 
tion, mixed regions of hot burned and cold unburned mate- 
rial emerge that exceed a cert ain length sc ale of the order 
cm toursi & Timm es 2006; Nie mever & WooslevI 
I1997I : iKhokhlov etaP 1 19971 : [Seitenzah l et al. 2009b^. Based on 
the Zel'dovich gradient mechanism (Zerdovich et aP I197Q|) a 
spontaneous ignition in such regions may lead to a suffi- 
ciently strong s hock for the formation of a detonatio n front 



jBlin n ikov & Khokhloy ||l9"8 ^ : lpokhlovlll991allb: Khokhl ov et sl\ 
I1997I : ISeitenzahl et all l2009al) . ILisewski et all OOOQO pointed out 
that v' has to excee d 10^ cms~^, but more recent studies of 
IWooslev et alj (l2009h suggested that smaller values of V ^ 5 x 
10^ cm s ^ may be sufficient to trigger the DDT. In three- 
dimensional simulations of pure deflagrations iRopk j bOOVh found 
that the probability P(V ^ 10^ cm s~^) of finding V of at least 
10^ cm s ^ may be high enough for a DDT to occur. 

The length scale 4rit is not resolved in our large-scale sim- 
ulations. Therefore a subgrid scale model is employed that eval- 
uates the probability P(V ^ 10^ cm s~^)(^) at the time t at spe- 
cific places at the deflagration flame that obey the following DDT 
constraints. To restrict our analysis to the immediate area of the 
flame, we define Xf^gi = 0.4-0.6 as the allowed fuel fraction (the 
carbon/oxygen composition) in the grid cells, which ensures that 
the level set propagates on average through the center of these 
cells. We further define pf^ei = (0.6-0.7) x 10^ g cm~^ as the al- 
lowed fuel densities in these cells, where t his interva l constitutes a 
subrange of the suggested fuel densities of lWooslevI 12007) where 
DDTs are expected to occur. The number N^^^^it) and length A(^) 
of the grid cells that meet the above-mentioned constraints define a 
flame surface area A^^^^^^ 

(0 - NL^Jt) X (A ( 0)^, where D is the 
frac tal dimen s ion. F ollowing iKersteinI (Il988h . ISreenivasanI (1199 ih 
and IWooslevI (l2007h we use D = 2.36. We define now A^Qi{t) = 
P{y' ^ 10^ cms"l)(0 X ^s part of the flame that is 

capable of performing a DDT. If Adet(0 exceeds a certain thresh- 
old 

^crit (we use Acrit = ^cr\v again with D = 2.36), it is checked 
whether this condition holds at least for half an eddy turnover 
time Tf 



eddy 



■ 0.5 X icrit/v' = 5 X 10~^ s. This is to ensure that 
fuel and ash become sufficiently mixed on the scale ^rit- If finally 
Adet(0 ^ ^crit holds for Teddy 1/2' detonations are initialized in the 
grid cells that contain the highest velocity fluctuations, where the 
number of ignitions is given by the current ratio Adet(0 ^^^crit- This 
DD T criterion is consistent with findin gs of microscopical stud- 
ies (IWooslevll2007l : IWooslev et alj|2009b and is largely resolution- 
independent (for further details see Ciaraldi-Schoolmann & Ropke, 
in preparation). 



2.4 Hydrodynamic evolution and explosion morphologies 

First we point out that all fourteen simulations result in a com- 
pletely gravitationally unbound remnant and are in this regard to be 
considered as "successful" explosions. The basic underlying evolu- 
tion of all models is rather similar to other extant three-dimensional 
delayed detonation simulations. 

T he burning fr ont propagates initially as a laminar conductive 
flame (iTimmes & Woos lev . l99Z) . The nuclear burning releases en- 
ergy by converting ^^C and ^^O to a mix of more tightly bound 
nuclei. This energy release is sufficient to at least partially lift the 
degeneracy of the burned material, which means that the ash is not 
only hotter but also less dense than the surrounding background 



material. Buoyancy forces then lead to a rise of the deflagration 
ash material against gravity. Once the deflagration bubble exceeds 
a critical size it becomes subjec t to the Rayleigh-T aylor (RT) insta- 
bility (e.g. Khokhlov 1995; To wnslev et alj[2007b , which wrinkles 
the flame front and causes a rapid acceleration of the growth of its 
surface area. As these RT-unstable plumes of hot deflagration ash 
continue to rise towards the stellar surface, the low Reynolds num- 
bers coupled with the rather high relative speeds lead to Kelvin- 
Helmholtz induced shear on the sides of the rising plumes. This 
shear-induced turbulence in turn results in the high velocity fluctua- 
tions that lead to the fulfillment of our DDT criterion (see Sec. 12. 31 ). 
The nascent detonation fronts propagate supersonically from their 
initiation sites and quickly burn any acc essible (i.e. topologically 
connected, cf . iMaier & Niemevej l2006h fuel that remains above 
the respective density threshold to IMEs and IGEs. As an aside, 
we mention that the detonation level set may be advected into low 
density (p < 5 x 10^ g cm~^) regions. However, this does not re- 
sult in any burning the re since we use th e fuel density dependent 
energy release tables of iFink et"aD (l201Qh . Once the detonation has 
burned all remaining and accessible fuel, the supernova enters the 
stage of ballistic expansion. In Fig.|2]we show the morphology for 
two representative models (N3 and NlOO) at three different stages 
of the evolution. 

Since we have essentially parametrized the strength of the 
deflagration phase by ignition kernel number (see Sec. 12.11 and 
Fig. [B, the models with the most ignition kernels burn the least 
amount of fuel in the detonation and vic e versa ( Ro pke et alj2007bl : 
iMazzaH et aDl2007l : iKasen et al.]l2009h . Furthermore, it is evident 
that models with fewer ignition kernels are more asymmetric. A 
prime example for this is N3, see left column of Fig.|2] The small 
number of ignition kernels necessarily leads to a highly asymmetric 
seed configuration of the RT-unstable deflagration plumes, which 
results in a highly asymmetric distribution of the deflagration ashes 
in the star. For our N3 case, only one hemisphere is burned by the 
deflagration. The fact that the detonation can only be triggered at 
the interface of the deflagration flame front and the fuel reinforces 
this asymmetry. The resulting distribution of the nucleosynthesis 
products fully reflects the asymmetric evolution of the deflagration 
and detonation (see Fig. [3]). In contrast, consider the NlOO model 
(right column in Fig.|2]). The rather symmetric ignition configura- 
tion results in deflagration plumes rising in all directions. As a con- 
sequence, DDTs also occur on all sides of the star and the explosion 
asymmetries are very modest, which again reflects in a more sym- 
metric distribution of burning products (see Fig.|4]). 



3 NUCLEOSYNTHESIS 

For each of the fourteen models, we have determined the iso- 
topic composition of the ejecta using ou r tracer particle metho d 
(iTravadio et al.ll2004l : iRopke et al.ll2006al : ISeitenzahl et"aDl201Qh . 
In each initial model we have distributed one million tracer par- 
ticles of equal mass in a way that the underlying density pro- 
file of the WD is reproduced by the tracer particles. Extrapolat- 
ing the tracer resolution and yield convergence study (done in 
2D) of ISeitenzahl et al.l fcOlOh to 3D, this number of tracer par- 
ticles (100 per axis) is sufficient to reliably predict the yields 
for the most abundant nuclides. For all models, the initial chem- 
ical composition for the post-processing is taken to be 47.5 per 
cent ^^C, 50 per cent ^^O, and 2.5 per cent ^^Ne by mass. The 
^^Ne fraction roughly corresponds to solar metallicity of the zero- 
age main- sequence progenitor under the often made approxima- 
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Figure 3. Tracer particle positions for the N3, N40, N300C, and N1600 models (from top to bottom). The tracer particles are colored by the mass fractions of 
^^Ni (first column), ^"^Fe+^^Fe+^^Ni (representing stable iron group nuclei, second column), ^^Si (third column), ^^O (forth column) and ^^C (fifth column) 
at r = 100 s. For ^^Ni, ^^Si, ^^O and the stable iron group nuclei, the tracer particle cloud is cut in the plane of the page and only the bottom hemisphere is 
shown. Plotting the data this way allows the viewer to also see the abundance distributions in the deep core while still retaining the three dimensional nature 
of the tracer particle locations. Only for ^^C the full spherical cloud is shown since ^^C largely resides near the surface. 



tion that ^^N is processed to ^^Ne during core helium burn- 
ing via i'^N(a,7)^^F(j3+)^^0(a,7)22Ne. The ^^Ne content intro- 
duces an excess of neutrons and results in an electron fraction 
of Yq = 0.49886. This composition is assumed to be homoge- 



an input to radiative transfer calculations (e.s 


[. Kasenet all 12006: 


Siml2007:IS 


m et al.l2010l:lKromer et alJl2010l 


IPakmor etalbOlO: 


Ropke et alJ 


2012h to derive model lisht curves and spectra. Fur- 



thermore, the yields can be used as an input for galactic chemical 
evolution calculations or to make predictions of the shapes of the 
late-time light curves re.g. lSeitenzahl et alJ2009d :[ Seitenzahll201ll : 
iRopke et al.|[2oT2h . 



3.1 Total integrated nucleosynthetic yields 

The yields of stable nuclides are presented in Table[2] For this table 
we have completely decayed all radioactive nuclides with half-lives 
less than 2 Gyr (such as ^^K or ^^Mn) to stability; nuclides with 
half-lives longer than 2 Gyr (i.e. ^^V) are tabulated with their pro- 
duction yields at ^ = 100 s. The yields at ? = 100 s of some abun- 
dant and long-lived radioactive nuclides are tabulated in Table [3] 
For convenience, we also show the yields of ^^Ni, IGEs, IMEs, 
1^0, and in Fig. [5] 

We note a few obvious trends in the total integrated yields. As 
expected, there is a clear trend that the total mass of ^^Ni produced 
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Figure 4. Tracer particle positions for the central density sequence of NlOO models (from the low central density NIOOL model in the top row to the high 
central density model NIOOH in the bottom row). Just as in Fig.[3]the tracer particles are colored by the mass fractions of ^^Ni (first column), ^^Fe+^^Fe+^^Ni 
(representing stable iron group nuclei, second column), ^^Si (third column), ^^O (forth column) and ^^C (fifth column) at f = 100 s. Again, for ^^C we show 
the whole cloud of tracer particles, whereas only one hemisphere is shown for all other species. 
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Figure 5. Fraction of the total mass at r = 100 s in ^^Ni, IGEs (includ- 
mg ^^Ni), IMEs, ^^O, and ^^C as a function of ignition kernel number. We 
use upward-pointing triangles for NIOOH, downward-pointing triangles for 
NIOOL, squares for N300C and N1600C, and circles for all other models. 



decreases with increasing ignition kernel number - from L 1 1 M© 
for Nl down to 0.32 M© for N1600C. T his covers the range o f 
expected ^^Ni masses of normal SNe la (IStritzinger et al 
Similarly, the total mass produced as intermediate mass elements 
(see e.g. ^^Si or ^^Si) shows a clear trend of increasing yields with 
increasing ignition kernel number. 

This can be understood in the following way. For small num- 
bers of ignition kernels, the deflagration is relatively weak and 
burns only a small fraction of the total mass of the star. The asso- 
ciated energy release results in only moderately strong expansion. 
Consequently, at the first occurrence of a DDT, the central density 
of the WD is still quite high and most of the remaining fuel will be 
burned to IGE (most of which is ^^Ni) by the detonation. In con- 
trast, the larger the number of ignition sparks, the stronger is the 
expansion during the deflagration phase. The low central density at 
the onset of the first DDT now results in much of the remaining fuel 
to be located below a density of 1 x 10^ g cm~^, approximately the 
cutoff where a detonation in equal mass carbon-oxygen material 
ceases to burn completely to NSE and instead produces IMEs as a 
result of incomplete burning. 

^^Ni, which is by far the most abundant stable nickel isotope in 
all models, shows remarkably little variation. The models produce 
^^Ni masses from 6.2 to 7.5 x 10~^ M©. The only exception is 
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the low central density model NIOOL, which undergoes less in situ 
neutronization and only synthesizes 3.8 x 10~^ of ^^Ni. 

While estimated ^^Ni masses of SNe la cover quite a range, 
most observed events cluster around 0.6 M© (e.g. Stritzin ger et al.1 
l2006h . For this reason, we have chosen the three NlOO models for 
a comparison of their isotopic iron peak production factors with 
the solar values, (^ Az/Xse^, ) / (- ^az/^5 6p. ) g) (see Fig. (6]). For the solar 
composition we use lLoddersI (l2003h . Our model production factors 
exhibit qualitatively the right behavior that is required for SN la 
yields: production factors of self-conjugate and slightly neutron 
rich isotopes ^^Cr, ^^Cr, ^^Cr, ^^Mn, ^^Fe, and ^^Ni only weakly 
depend on the particular choice of central density of the WD at the 
time of ignition. These isotopes are expected to be largely synthe- 
sized in nuclear statistical equilibrium in SNe la, which requires 
production factors > 1 as observed here. Large overproduction fac- 
tors ([^z/56pg] > 2)^ are marginally in conflict with the requirement 
that S Ne la produced roughly half of the ^^Fe present in the Sun 
today (IClavtonll2003h . We note, however, that models with differ- 
ent central density at ignition yield lower overproduction factors 
(Fig. (6]) and that lower progenitor metallicities also results in less 
neutron rich iron group isotopes (see Section [T2l ). 

The only isotope with an overproduction factor > 3 is ^^Cr for 
the high central density model NIOOH. We interpret the large over- 
production factor we obtain there as an indication that delayed det- 
onation SNe la that ignite at central densities > 5.5 x 10^ g cm~^ 
are rare and constitute at most a small fraction of all SN la events. 
The most neutron-rich stable Fe-peak isotopes ^^Cr, ^^Fe, and ^^Ni 
are shielded by ^^Fe, ^^Ni, and ^^Zn from thQ Z = N line and 
thus require the most neutronization for direct production. Con- 
sequently, these isotopes are the most sensitive to the central den- 
sity, with the highest production factors occurring for the NIOOH 
model. The very pronounced underproduction of ^^Fe and ^^Ni in 
all the models is not a problem, since the S-process is the dominant 
source of nucleosynthesis for these isotopes. The solar abundances 
of the remaining isotopes ^^Co, ^^Ni, ^^Ni, and ^^Ni contain signif- 
icant contributions from explosive Si-burning (core collapse SNe), 
alpha-rich freeze-out (including SNe la), and the S-process. Our 
production factors for SNe la on the order of a few to several tens 
of percent are therefore also very reasonable. 

For a detailed discussion of the nucleosynth esis origi n in th e 
Sun of the isotopes discussed here see the book bv lClavtonl ( [2003h . 
which we have used as a reference for our comparison. 

The subset of our three NlOO models contains another inter- 
esting trend. We find that the ^^Ni mass in the NlOO model se- 
quence increases with central density, i.e. the high central density 
model NIOOH produces the most and the low central density model 
NIOOL produces the lea st amount of ^^Ni. Th is is exactly opposite 
to the behavior found bv lKruegeretaPfeOld) when they analyzed 
their statistical sample of 150 two-dimensional delayed-detonation 
explosion simulations. They considered a range of central ignition 
densities (1 to 5 x 10^ g cm~^), a fixed deflagration to detona- 
tion transition density of 10^-^ g cm~^, but did not perform post- 
processing with a detailed nuclear reaction network. 

In our models, the fracti on of NSE material (for a recent dis- 
cussion of NSE in SN la see ISeitenzahl et al.|[2009dh that is pro- 
duced as ^^Ni is highest for the low central density model and 
lowest for the high central density model, which is what is ex- 
pected. The total amount of iron group material synthesized is, 
however, strongly incr easing with central dens ity. We have seen the 
same trend already in lSeitenzahl et al.l (l201lh . where we had esti- 
mated a roughly constant ^^Ni yield as a function of central den- 
sity. The ^^Ni mass in that work was derived from the total mass 
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Figure 6. Production factors (^^z/x^^^J / {Xaz/x5(,p^)q of stable iron peak 
isotopes relativ e to ^^Fe normalized to the corresponding solar values 
(lLoddersll2003l) for models NIOOL (downward-pointing triangles/dashed 
lines), NlOO (circles/solid lines), and NIOOH (upward-po inting trian- 
gles/do tted lines) and W7 (transparent squares/thin solid lines i lMaeda et al.l 
bOlObh . The lower the central density, the smaller the overproduction of 
^^Fe and ^^Ni. 

in IGE and the electro n fraction Yq, essentially using the formula 
in lTimmes et al.1 (l2003h . Here we determine the nucleosynthesis in 
detail with one million tracer particles and we find that the ^^Ni 
mass increases with increasing central density. We caution, how- 
ever, that our sample consisting of three models in a single ignition 
configuration is not statistically significant. 

3.2 Dependence of yields on progenitor metallicity 

To assess the impact of varying the progenitor metallicity, we have 
also post-processed the NlOO model with one-half, one tenth, and 
one hundredth of the canonical ^^Ne mass fraction of 0.025. The 
^^C mass fractions were thus 0.4875, 0.4975 and 0.49975 and 
was kept constant at 0.5. The models with reduced ^^Ne are 
called N100_Z0.5, N100_Z0.1, and N100_Z0.01 respectively and 

their yields are als o presented in Tables [ 3 and [3] 

As expected (iTimmes et al.ll2003l : Fl7avaglio et al.ll2005h . the 
^^Ni yields increases with decreasing initial ^^Ne due to the de- 
creasing electron fraction Fe, largely at the cost of stable iron group 
isotopes such as ^^Ni or ^^Fe. As a consequence, the overproduc- 
tion factors of these isotopes are successively reduced to < 2 even 
for the canonical central density case (Fig. IT}. Since the ^^Fe in 
the Sun is largely due to supernovae that had progenitors with sub- 
solar metallicity, we argue that the isotopic Fe-group yields of our 
delayed-detonations are not inconsistent with solar isotopic ratios. 

3.3 Spatial and velocity distribution of nucleosynthetic yields 

Total integrated yields constitute very important predictions and di- 
agnostics for an explosion model. How the yields are distributed in 
mass and velocity space is, however, at least equally important, es- 
pecially when it comes to direct observables of the supernova such 
as spectra or light curves. 

It is not possible to show the spatial yield distributions of all 
nuclei for all models. Therefore, we restrict ourselves to show only 
^^Ni (as the main source of radioactive heating), ^^Fe+^^Fe+^^Ni 
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Figure 7. Production factors I {^^z/^5(>yiq)q of stable iron 

peak isotope s relativ e to ^ ^Fe normalized to the corresponding so- 
lar values ( iLoddersI l2003h for models NlOO (circles/solid lines), 
N100_Z0.5 (diamonds/dot-dashed lines), N100_Z0.1 (upward-pointing tri- 
angles/dashed lines), N100_Z0.01 (downward-poi nting triangles/dotted 
lines) and W7 (transparent squares/thin solid lines: iMaeda et al.ll2010bh . 
The smaller the initial ^^Ne mass fraction, the smaller the overproduction 

(as most abundant stable IGE s and important sources of opacity), 
^^Si (as the generally most abundant intermediate mass isotope) 
and the fuel isotopes ^^O and ^^C, for a representative sample of 
the four models N3, N40, N300C, and N1600 (see Fig.©. Globally, 
the abundances are stratified as is expected for delayed-detonations: 
Near the surface some unburned carbon sits on top of an oxygen- 
rich layer that is composed of unburned oxygen fuel and products 
of low-density carbon burning. Further inwards IMEs such as ^^Si 
and ^^S are the most abundant species, although a strict onion- shell 
structure as is seen in one-dimensional models is not present. In- 
stead, Rayleigh-Taylor plumes of material rich in stable IGEs have 
penetrated into this layer. Consequently, stable IGEs are not found 
at th e lowest velocities as predicted by one-dimensional model s 
(e.g. iNomoto et all 1 19841 : iKhokhlovl Il99ld : iHoflich et alJ l2004h . 
but rather at intermediate velociti es (^3,000— 10,000 kms~^). 
This agrees well with the results of iMaeda et al.1 (l2010bh . who find 
typical velocities of 5,000 — 10,000 km s~^ for the deflagration 
ash in their two-dimensional 0-DDT model. This signature of the 
burning in the deflagration is characteristic to multi-dimensional 
simulations that are not ignited at the very center. Unlike in the 
one-dimensional case, the hot, less dense, buoyant ash can float 
towards the surface. We note that Stehle et al. (2005) find substan- 
tial amounts of stable iron out to velocities of about 9,000 analyz- 
ing spectra of SN 2002bo. In contrast to our model of comparable 
brightness (NlOO), in their abundance tomography, which is based 
on one-dimensional explosion models, the abundance of stable iron 
increases steadily towards lower velocities and dominates in the 
very center. 

In between the IGE-rich plumes, small pockets filled with 
oxygen-rich material remain - these were downdrafts in the defla- 
gration that were burned only in the latest phases of the detonation 
at low densities. The central regions, on the other hand, which were 
burned by the detonation at high densities to NSE, form a homo- 
geneous ^^Ni clump. In the fainter models, this clump fragments 
progressively. This is a natural consequence of the fact that in these 
models the majority of IGE material is produced by the deflagra- 



tion and not by the detonation. In these models, the strong initial 
expansion due to the deflagration results in such low densities in 
the core that the ensuing detonation fails to process this material 
to NSE (see Fig. [3jp)). We also show the same set of mass frac- 
tions for the three NlOO models (see Fig.|4l). The global picture of 
chemical stratification discussed above is rather insensitive to cen- 
tral ignition density. 

To better visualize the underlying main trends of the whole 
model suite, and to facilitate a comparison of the three-dimensional 
models, we reduce the information and show one-dimensional 
abundance profiles in velocity space (see Fig. [8]). We emphasize 
that this averaging or binning of the three dimensional data erases 
all information about the inhomogeneities and (sometimes pro- 
nounced) asymmetries brought about by e.g. the rising plumes of 
deflagration ash or downdrafts of nuclear fuel. This is also true 
for models with a strong (turbulent) deflagration phase, which are 
rather symmetric under rotations on large scales, but exhibit strong 
inhomogeneities in the burning products on small scales (see the 
different morphology plots in Sec. 12.41 ). For example, the presence 
of e.g. ^^O and ^^Ni in these one-dimensional profiles at the same 
velocity does not necessarily imply a co-spatial existence of these 
nucl ear species. 

iKozma et"aD (l2005h showed that ejecta with oxygen at low 
velocities (as we find it for our models with large A^^) may lead 
to strong [Ol] AA6300,6364 emission at late times, which is not 
observed in SNe la. The mere presence of oxygen at low veloci- 
ties, however, is not the only relevant condition for [O I] emission. 
Whether [O I] features will arise, depends strongly on the ionisa- 
tion state and a possible microscopic m ixing of different sp ecies. 
Compared to the deflagration model of iKozma et al.l (l2005h . our 
delayed detonation models have lower density ejecta, which could 
lead to a higher ionization. The question of microscopic mixing 
cannot be answered from present-day numerical models since it 
takes place on scales which are not resolved. However, if such a 
mixing is present, stronger transitions than [O I] will dominate the 
cooling. 

As a general trend, ^^Ni is hardly present at velocities above 
~ 12, 000 kms~^Q The exceptions are the models N3 and N5, 
where parts of the asymmetrically rising deflagration plumes have 
already risen to the stellar surface when the first DDT occurs (see 
Fig Efc)- In these models, the deflagration ash is therefore not com- 
pletely enclosed by burning products of the detonation and there is 
some ^^Ni present at very high velocity. ^^Ni and ^^Ni are produced 
co-spatially with ^^Ni and more or less follow its distribution, albeit 
with a lower abundance. 

The distribution of IMEs in velocity space shows a clear trend 
with ignition kernel number. The brightest models (small number 
of ignition kernels) contain IMEs from the highest expansion veloc- 
ities down to ~6, 000-7, 000 km s~^ . For the fainter models (larger 
number of ignition kernels) the inner boundary of the IMEs contin- 
ues to move inward and the IMEs at the highe st velocities are mor e 
and more replaced by unburned fuel (see also lMazzali et al.ll2007h . 

Notably, carbon and oxygen are extending down to low ve- 
locities. There is always some carbon present down to or even be- 
low velocities of 12, 000 km s ~^ . Note that this is much lower than 
what has been found for W7 jNomoto et al.lll984l : llwamoto et alj 
ll999l : lMaeda et al. 11201 Obh . Oxygen reaches down even farther, al- 
ways present below ^8,000 km s~^ and sometimes even reaching 



One should be aware of the fact that already traces of IGEs on the equiv- 
alent level of solar abundance may affect the observables significantly. 
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to the very center. The presence of unburned fuel in our models is a 
distinct multi-dimensional effect. For fainter models, the products 
of low-density detonation burning extend down to lower velocities. 
In extreme cases, pockets of unburned fuel remain near the center 
that are not reached by the detonation at all. 

The other abundant iron group isotopes ^^Co, ^^Fe, and ^^Fe 
are mainly produced in the deflagration phase. As a result, they 
are concentrated in an off center shell surrounding the central ^^Ni 
bubble (see Figs.[3}|8l). Only models with a weak deflagration phase 
(in particular Nl, N3, N5) exhibit high enough central densities 
during the detonation phase that some neutronization and normal 
freeze-out from NSE also synthesizes significant amounts of these 
isotopes in a second, centrally concentrated production site (see 
Fig. [8}. 

The idea suggests itself to compare the yield morphologies of 
our three-dimensional simulations to existing two-dimensional re- 
sults from the literature. Such a comparison, ho wever, is difficult to 
make quantitatively. For example, the results of iKasen et al ]([200i) 
were not based on detailed yield information obtained by post- 
processing each model but mapped representative yields into the 
individual realizations. A qualitative discussion, however, is possi- 
ble: Comparing Fig.[3]to their figure 1, we can, however, identify a 
key commonality. In both sets, the stable IGEs in bright models are 
most abundant at intermediate velocities, surrounding a core rich 
in ^^Ni. Conversely, the fainter models have the stable IGE much 
closer to the center. To the best of our knowledge, the only pub- 
lished results of multi-dimensional delayed-detonati ons with de- 
tailed n ucleosynthetic yields are the two models of iMaeda et al.l 
boiobh . A meaningful comparison to their C-DDT model, which 
is based on a rather symmetrically and centrally ignited deflagra- 
tion, is not possible. In their C-DDT model, the DDT occurs at 
such a late time and low density that no more ^^Ni is produced by 
the detonation wave, resulting in a very low mass of final ^^Ni of 
~O.25M0. This is both outside the range of^^Ni masses covered by 
our models as well as observed "normal" SNe la. On the other hand, 
their 0-DDT model, which is based on an asymmetric ignition of 
the deflagration, produces ~O.54M0 of ^^Ni, which is rather sim- 
ilar to our N150 or NlOO models. In spite of the different symme- 
tries in the ignition configuration, the angle average d distribution 
of nucli des in velocity space (compare figure 11(c) of iMaeda et all 
bOlObh and Fig.H are quite similar. In both 0-DDT and our mod- 
els, the lowest velocities are dominated by ^^Ni and ^^Ni, IMEs are 
present down to velocities of a few thousand km s~^, ^^O down 
to just below 5,000 km s~^, and unburned ^^C down to around 
^10, 000 km s~^ . Stable IGE are only found at intermediate veloc- 
ities between ~5,000 - 11,000 km s~^ and ^^Ni cuts off sharply 
around 11,000 km s-^ 

In 2D axisymmtric simulations, the rising deflagration "bub- 
bles" are in fact "tori". This leads to differences in the morphology, 
even if angle averages look alike. 3D models allow for variation on 
smaller angular scales. A rising torus is a more symmetric and also 
more extended object than a rising bubble. Since it contains much 
more volume than a single bubble, a 2D deflagration starting from 
the same ignition site releases more energy and is more symmetric 
than the corresponding case in 3D. This may explain why the 2D 
O-DDT model, which was ignited in a rather asymmetric configu- 
ration in only 29 spots, ends up looking rather symmetric and not 
nearly as bright as e.g. our N20 or N40 model. 



4 SUMMARY AND CONCLUSIONS 

We have performed fourteen three-dimensional hydrodynamical 
simulations for delayed detonation SNe la for a range of ignition 
conditions. For each simulation we have determined the complete 
nucleosynthetic yields by post-processing one million tracer par- 
ticles with a nuclear reaction network. This set of models consti- 
tutes the first suite of three-dimensional explosion models that cov- 
ers the range of expected ^^Ni masses of spectroscopically normal 
SNe la. From our study we conclude that, fixing all other parame- 
ters of the exploding WD but the ignition configuration, a delayed- 
detonation of "normal" SN la brightness likely requires rather sym- 
metrical, central ignition to occur in nature. Only such a setup suf- 
ficiently pre-expands the WD in the deflagration phase to reduce 
the ^^Ni production in the subsequent detonation. Otherwise only 
the brightest SNe la could be explained with delayed detonations 
of Chandrasekhar-mass WDs as realized in our models by few (and 
thus asymmetrically distributed) ignition sparks. W hether or not 
the dem and for symmetric ignition is in conflict with lMaeda et alj 
(l2010ah , who require a typical off-set of ~3, 500 km s~^ in the de- 
flagration ashes, is not clear. In spite of rather symmetric ignition 
in our "normal" models, the fact that one or a few RT-modes grow 
faster than the others still results in a distribution of stable Fe-peak 
isotopes that is not uniform (see second column of Fig.|4]). 

Our suite of models is the first published set of three- 
dimensional delayed detonation simulations with detailed isotopic 
nucleosynthetic yields. As such, the yields presented here lend 
themselves to be used as an input for Galactic chemical evolution 
calculations. They also set the stage for predicting observables by 
radiative transfer calculations. Overall, we expect the brightness 
range of normal SNe la to be covered by our set of models. It re- 
mains to be seen whether the predicted spectra match the observa- 
tions and whether the set of models follows the width-luminosity 
relation (iPhillip s 1993; Ph illips et al. 1999 ) and other observational 
trends, as a set of two-dimensional delay ed-detonation models did 
(iKasen et all l2009l : iBlondin et all 120111) . The Fe-group isotopes 
in our models of normal SN la brightness (^^Ni masses around 
0.6 Mq) are synthesized in the required proportions (see Figs.[6]and 
|7]), which tells us that delayed-detonations cannot be ruled out as 
the dominant SN la explosion channel based on solar isotopic Fe- 
group ratios. The presence of IGEs (in particular stable isotopes) 
at high velocities and oxygen and carbon in the inner ejecta in our 
models is expected to leave testable imprints on the observables. 
Moreover, our set of three-dimensional full- star models captures 
asymmetries that potentially could be constrai ned by spectrapo- 
larim etry measurements (see e.g. the review bv lWang& Wheeled 
.2008.) . 
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velocity [10^ cm s'^ ] 



Figure 8. Shown are average mass fractions of some select isotopes and IME at time t ^ 100 s in velocity space. 
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Table 2. Asymptotic nucleosynthetic yields (in solar masses) of stable nuclides. 





Nl 


N3 


N5 


NIO 


N20 


N40 


NIOOH 


NlOO 


NIOOL 


N150 


N200 


N300C 


N1600 


N1600C 


N100_Z0.5 


N100_Z0.1 


N100_Z0.01 




2.61E-03 


9.90E-03 


9.05E-03 


4.43E-03 


9.20E-03 


3.90E-03 


3.87E-03 


3.04E-03 


3.85E-03 


1.72E-02 


1.21E-02 


8.86E-03 


1.06E-02 


1.68E-02 


3.10e-03 


3.15e-03 


3.16e-03 




1.84E-08 


6.15E-08 


5.05E-08 


2.57E-08 


4.52E-08 


2.18E-08 


2.28E-08 


1.74E-08 


2.17E-08 


l.OOE-07 


6.57E-08 


4.95E-08 


5.57E-08 


8.44E-08 


8.47e-09 


1.91e-09 


2.72e-10 


14n 


2.92E-06 


9.93E-06 


8.46E-06 


3.85E-06 


8.34E-06 


4.30E-06 


4.25E-06 


3.21E-06 


3.98E-06 


1.84E-05 


1.33E-05 


9.16E-06 


1.04E-05 


1.88E-05 


1.80e-06 


4.71e-07 


7.22e-08 


15n 


3.36E-09 


1.22E-08 


1.03E-08 


4.47E-09 


1.03E-08 


5.16E-09 


5.24E-09 


3.67E-09 


4.66E-09 


2.29E-08 


1.71E-08 


1.14E-08 


1.29E-08 


2.41E-08 


2.07e-09 


2.98e-09 


8.73e-08 


16o 


2.63E-02 


4.74E-02 


5.63E-02 


5.16E-02 


9.04E-02 


9.89E-02 


7.30E-02 


l.OlE-01 


1.24E-01 


1.24E-01 


1.96E-01 


1.21E-01 


1.91E-01 


2.72E-01 


9.87e-02 


9.64e-02 


9.47e-02 


1^0 


3.96E-07 


1.37E-06 


1.16E-06 


5.34E-07 


1.12E-06 


5.54E-07 


5.61E-07 


4.13E-07 


5.14E-07 


2.48E-06 


1.74E-06 


1.22E-06 


1.36E-06 


2.42E-06 


2.84e-07 


9.32e-08 


5.43e-09 


18o 


3.32E-09 


1.33E-08 


l.llE-08 


4.54E-09 


1.12E-08 


5.29E-09 


5.59E-09 


3.53E-09 


4.61E-09 


2.52E-08 


1.98E-08 


1.25E-08 


1.44E-08 


2.73E-08 


2.23e-09 


1.21e-09 


9.93e-10 


19f 


3.73E-11 


1.35E-10 


1.17E-10 


5.05E-11 


1.22E-10 


6.22E-11 


6.32E-11 


4.39E-11 


5.68E-11 


2.64E-10 


2.13E-10 


1.36E-10 


1.58E-10 


2.97E-10 


2.20e-ll 


1.48e-ll 


4.79e-ll 


20Ne 


1.47E-03 


3.37E-03 


3.75E-03 


2.40E-03 


5.41E-03 


4.15E-03 


3.66E-03 


3.53E-03 


4.33E-03 


8.72E-03 


1.15E-02 


6.76E-03 


9.40E-03 


1.73E-02 


3.60e-03 


3.69e-03 


3.74e-03 


21 Ne 


3.08E-07 


9.79E-07 


8.81E-07 


4.16E-07 


9.63E-07 


5.43E-07 


5.20E-07 


4.11E-07 


5.17E-07 


1.98E-06 


1.68E-06 


1.08E-06 


1.29E-06 


2.39E-06 


1.97e-07 


4.47e-08 


6.93e-09 



22 Ne 


6.40E-05 


3.26E-04 


2.87E-04 


1.31E-04 


2.58E-04 


5.77E-05 


7.62E-05 


4.07E-05 


5.51E-05 


4.83E-04 


2.34E-04 


2.11E-04 


2.32E-04 


2.97E-04 


1.65e-05 


2.30e-06 


1.71e-07 


23 Na 


2.20E-05 


6.41E-05 


6.09E-05 


3.22E-05 


7.30E-05 


4.66E-05 


4.25E-05 


3.74E-05 


4.68E-05 


1.38E-04 


1.38E-04 


8.53E-05 


1.09E-04 


2.02E-04 


2.63e-05 


1.96e-05 


1.72e-05 


24Mg 


3.93E-03 


7.13E-03 


8.53E-03 


7.77E-03 


1.46E-02 


1.54E-02 


1.15E-02 


1.52E-02 


1.83E-02 


1.93E-02 


3.32E-02 


1.97E-02 


3.08E-02 


4.64E-02 


2.02e-02 


2.69e-02 


2.90e-02 


25 Mg 


3.35E-05 


9.26E-05 


8.92E-05 


5.07E-05 


l.llE-04 


7.70E-05 


6.86E-05 


6.49E-05 


8.02E-05 


2.02E-04 


2.14E-04 


1.33E-04 


1.75E-04 


3.12E-04 


3.09e-05 


1.06e-05 


8.99e-07 


26 Mg 


5.15E-05 


1.36E-04 


1.34E-04 


7.55E-05 


1.71E-04 


1.17E-04 


1.04E-04 


9.66E-05 


1.19E-04 


3.07E-04 


3.27E-04 


2.01E-04 


2.63E-04 


4.82E-04 


4.44e-05 


7.36e-06 


1.04e-06 



27 Al 


1.98E-04 


3.95E-04 


4.56E-04 


3.71E-04 


7.32E-04 


7.05E-04 


5.47E-04 


6.74E-04 


8.32E-04 


1.04E-03 


1.64E-03 


9.68E-04 


1.48E-03 


2.37E-03 


5.88e-04 


2.68e-04 


8.71e-05 


28si 


6.32E-02 


8.99E-02 


1.19E-01 


1.38E-01 


1.98E-01 


2.59E-01 


2.12E-01 


2.84E-01 


3.55E-01 


2.71E-01 


3.28E-01 


3.19E-01 


3.61E-01 


3.44E-01 


2.90e-01 


2.94e-01 


2.89e-01 


29Si 


2.69E-04 


4.64E-04 


5.68E-04 


5.17E-04 


9.49E-04 


1.03E-03 


7.73E-04 


1.03E-03 


1.25E-03 


1.30E-03 


2.08E-03 


1.29E-03 


1.97E-03 


2.86E-03 


7.28e-04 


4.30e-04 


1.35e-04 


30Si 


5.96E-04 


l.OOE-03 


1.23E-03 


1.18E-03 


2.12E-03 


2.35E-03 


1.72E-03 


2.36E-03 


2.86E-03 


2.77E-03 


4.76E-03 


2.87E-03 


4.53E-03 


6.55E-03 


1.19e-03 


1.44e-04 


1.84e-05 


31p 


1.40E-04 


2.28E-04 


2.87E-04 


2.78E-04 


4.87E-04 


5.60E-04 


4.20E-04 


5.77E-04 


7.09E-04 


6.59E-04 


1.08E-03 


6.85E-04 


1.05E-03 


1.47E-03 


3.58e-04 


1.05e-04 


3.54e-05 



32s 


2.62E-02 


3.70E-02 


4.79E-02 


5.74E-02 


7.74E-02 


l.OlE-01 


8.55E-02 


l.llE-01 


1.38E-01 


1.07E-01 


l.lOE-01 


1.27E-01 


1.22E-01 


1.03E-01 


1.12e-01 


1.12e-01 


1.15e-01 




7.51E-05 


1.06E-04 


1.42E-04 


1.53E-04 


2.43E-04 


3.14E-04 


2.37E-04 


3.39E-04 


4.21E-04 


3.27E-04 


5.23E-04 


3.65E-04 


5.47E-04 


6.79E-04 


2.39e-04 


1.04e-04 


4.57e-05 


34s 


8.26E-04 


1.16E-03 


1.57E-03 


1.75E-03 


2.84E-03 


3.73E-03 


2.74E-03 


4.04E-03 


5.02E-03 


3.68E-03 


6.22E-03 


4.22E-03 


6.56E-03 


8.06E-03 


1.86e-03 


2.60e-04 


7.14e-06 


36s 


8.12E-08 


1.35E-07 


1.65E-07 


1.41E-07 


2.54E-07 


2.57E-07 


1.89E-07 


2.47E-07 


3.05E-07 


3.59E-07 


5.42E-07 


3.23E-07 


4.97E-07 


7.68E-07 


3.86e-08 


1.64e-09 


1.73e-ll 


35C1 


4.29E-05 


6.20E-05 


8.27E-05 


8.37E-05 


1.35E-04 


1.67E-04 


1.27E-04 


1.78E-04 


2.27E-04 


1.89E-04 


2.89E-04 


2.00E-04 


2.98E-04 


3.84E-04 


9.91e-05 


2.64e-05 


5.65e-06 



37 CI 


7.32E-06 


9.62E-06 


1.34E-05 


1.53E-05 


2.26E-05 


3.14E-05 


2.44E-05 


3.51E-05 


4.49E-05 


3.14E-05 


4.66E-05 


3.63E-05 


5.22E-05 


5.75E-05 


2.27e-05 


9.14e-06 


3.52e-06 


36 Ar 


4.52E-03 


6.36E-03 


8.03E-03 


9.89E-03 


1.28E-02 


1.61E-02 


1.43E-02 


1.77E-02 


2.17E-02 


1.76E-02 


1.50E-02 


2.08E-02 


1.64E-02 


1.23E-02 


1.85e-02 


1.92e-02 


2.04e-02 


38 Ar 


3.77E-04 


5.15E-04 


7.18E-04 


8.01E-04 


1.25E-03 


1.72E-03 


1.29E-03 


1.91E-03 


2.48E-03 


1.69E-03 


2.69E-03 


1.96E-03 


2.98E-03 


3.42E-03 


8.31e-04 


1.19e-04 


5.40e-06 


40 Ar 


1.68E-09 


2.45E-09 


3.13E-09 


2.90E-09 


4.79E-09 


5.27E-09 


3.81E-09 


5.21E-09 


6.60E-09 


6.64E-09 


1.04E-08 


6.29E-09 


9.87E-09 


1.48E-08 


4.66e-10 


9.87e-12 


5.11e-14 


39k 


2.26E-05 


2.98E-05 


4.12E-05 


4.65E-05 


6.80E-05 


9.52E-05 


7.40E-05 


1.07E-04 


1.39E-04 


9.53E-05 


1.42E-04 


l.lOE-04 


1.60E-04 


1.75E-04 


6.25e-05 


1.89e-05 


3.57e-06 



41k 


1.29E-06 


1.65E-06 


2.30E-06 


2.66E-06 


3.80E-06 


5.38E-06 


4.22E-06 


6.08E-06 


7.87E-06 


5.36E-06 


7.81E-06 


6.21E-06 


8.93E-06 


9.65E-06 


3.85e-06 


1.42e-06 


4.92e-07 


40Ca 


4.05E-03 


5.74E-03 


7.09E-03 


8.82E-03 


1.13E-02 


1.35E-02 


1.24E-02 


1.47E-02 


1.75E-02 


1.50E-02 


1.07E-02 


1.78E-02 


l.lOE-02 


7.50E-03 


1.57e-02 


1.66e-02 


1.77e-02 


42Ca 


8.32E-06 


l.llE-05 


1.55E-05 


1.74E-05 


2.65E-05 


3.69E-05 


2.80E-05 


4.15E-05 


5.44E-05 


3.63E-05 


5.71E-05 


4.22E-05 


6.44E-05 


7.30E-05 


1.67e-05 


2.08e-06 


8.87e-08 


43Ca 


1.63E-08 


1.93E-08 


2.63E-08 


2.92E-08 


4.29E-08 


5.62E-08 


4.37E-08 


6.12E-08 


7.44E-08 


5.57E-08 


8.61E-08 


6.03E-08 


9.21E-08 


1.14E-07 


2.35e-08 


6.22e-09 


6.23e-09 


44Ca 


3.25E-06 


4.40E-06 


5.24E-06 


6.49E-06 


8.14E-06 


9.14E-06 


8.92E-06 


l.OOE-05 


l.llE-05 


1.04E-05 


6.45E-06 


1.19E-05 


6.08E-06 


4.06E-06 


1.08e-05 


1.18e-05 


1.27e-05 



46ca 


1.58E-11 


1.99E-11 


2.66E-11 


2.64E-11 


4.01E-11 


4.76E-11 


7.13E-10 


4.85E-11 


6.06E-11 


5.50E-11 


8.60E-11 


5.65E-11 


8.43E-11 


1.24E-10 


2.24e-12 


5.88e-13 


5.82e-13 


48Ca 


1.97E-15 


3.88E-15 


7.61E-15 


7.56E-15 


1.86E-14 


5.84E-15 


9.97E-11 


7.85E-15 


3.04E-15 


6.17E-15 


5.21E-15 


8.89E-14 


4.52E-15 


6.99E-15 


5.04e-15 


4.85e-15 


4.82e-15 


45Sc 


5.12E-08 


6.53E-08 


8.56E-08 


1.03E-07 


1.36E-07 


1.83E-07 


1.56E-07 


2.05E-07 


2.51E-07 


1.92E-07 


2.20E-07 


2.20E-07 


2.48E-07 


2.45E-07 


1.43e-07 


7.69e-08 


3.94e-08 


46Ti 


4.33E-06 


5.72E-06 


8.01E-06 


9.13E-06 


1.33E-05 


1.87E-05 


1.46E-05 


2.11E-05 


2.75E-05 


1.86E-05 


2.72E-05 


2.17E-05 


3.11E-05 


3.34E-05 


9.34e-06 


1.37e-06 


2.19e-07 


47Ti 


1.80E-07 


2.28E-07 


3.01E-07 


3.58E-07 


4.84E-07 


6.33E-07 


5.44E-07 


7.05E-07 


8.51E-07 


6.53E-07 


7.66E-07 


7.44E-07 


8.46E-07 


8.59E-07 


3.37e-07 


9.30e-08 


9.90e-08 



48Ti 


1.02E-04 


1.41E-04 


1.67E-04 


2.07E-04 


2.67E-04 


2.85E-04 


2.87E-04 


3.14E-04 


3.36E-04 


3.31E-04 


1.80E-04 


3.78E-04 


1.50E-04 


9.30E-05 


3.43e-04 


3.73e-04 


3.87e-04 


49^1 


8.06E-06 


1.07E-05 


1.29E-05 


1.61E-05 


2.07E-05 


2.35E-05 


2.35E-05 


2.59E-05 


2.76E-05 


2.69E-05 


1.72E-05 


3.05E-05 


1.62E-05 


1.03E-05 


2.06e-05 


1.27e-05 


6.83e-06 


50Ti 


2.48E-08 


5.52E-08 


1.33E-07 


1.44E-07 


3.43E-07 


1.62E-07 


1.02E-04 


2.60E-07 


6.88E-10 


1.61E-07 


4.58E-08 


1.25E-06 


1.19E-09 


1.72E-09 


2.52e-07 


2.47e-07 


2.46e-07 


SOy 


8.95E-10 


1.21E-09 


1.94E-09 


2.34E-09 


3.70E-09 


4.68E-09 


3.20E-08 


5.56E-09 


3.84E-09 


5.07E-09 


6.16E-09 


6.72E-09 


6.11E-09 


7.52E-09 


3.01e-09 


2.30e-09 


2.28e-09 


51v 


3.41E-05 


4.14E-05 


5.01E-05 


6.18E-05 


8.06E-05 


9.51E-05 


1.51E-04 


1.05E-04 


l.OlE-04 


1.09E-04 


8.19E-05 


1.20E-04 


8.09E-05 


6.01E-05 


7.65e-05 


3.99e-05 


2.62e-05 



50Cr 


l.OlE-04 


1.39E-04 


1.84E-04 


2.27E-04 


3.08E-04 


4.06E-04 


3.95E-04 


4.51E-04 


4.63E-04 


4.47E-04 


4.48E-04 


5.05E-04 


5.08E-04 


4.01E-04 


2.53e-04 


l.lle-04 


8.59e-05 


52Cr 


3.55E-03 


4.06E-03 


4.56E-03 


5.45E-03 


7.23E-03 


7.70E-03 


1.16E-02 


8.57E-03 


8.02E-03 


8.90E-03 


5.62E-03 


9.85E-03 


4.78E-03 


3.94E-03 


8.97e-03 


9.44e-03 


9.64e-03 


53Cr 


5.25E-04 


5.27E-04 


5.91E-04 


7.09E-04 


9.48E-04 


1.12E-03 


1.46E-03 


1.24E-03 


9.98E-04 


1.29E-03 


1.06E-03 


1.36E-03 


1.06E-03 


9.20E-04 


1.02e-03 


7.61e-04 


6.27e-04 


54Cr 


9.85E-07 


2.06E-06 


4.54E-06 


5.18E-06 


1.08E-05 


7.59E-06 


1.25E-03 


2.61E-07 


l.llE-05 


8.55E-06 


4.94E-06 


3.24E-05 


3.17E-06 


4.06E-06 


1.07e-05 


1.05e-05 


1.05e-05 


55 Mn 


8.74E-03 


6.93E-03 


6.80E-03 


7.61E-03 


9.70E-03 


1.21E-02 


1.39E-02 


9.29E-03 


1.33E-02 


1.39E-02 


1.31E-02 


1.35E-02 


1.31E-02 


1.13E-02 


l.lle-02 


8.73e-03 


7.84e-03 



54Fe 


3.02E-02 


3.07E-02 


3.84E-02 


4.77E-02 


6.71E-02 


9.00E-02 


9.80E-02 


9.94E-02 


5.70E-02 


1.04E-01 


1.13E-01 


1.03E-01 


1.22E-01 


l.OOE-01 


8.19e-02 


6.87e-02 


6.62e-02 


56Fe 


l.llE-00 


1.04E-00 


9.80E-01 


9.46E-01 


7.90E-01 


6.70E-01 


7.47E-01 


6.22E-01 


5.33E-01 


5.84E-01 


4.36E-01 


5.33E-01 


3.90E-01 


3.54E-01 


6.46e-01 


6.66e-01 


6.72e-01 


57Fe 


3.37E-02 


3.05E-02 


2.91E-02 


2.82E-02 


2.31E-02 


2.03E-02 


2.25E-02 


1.88E-02 


1.43E-02 


1.79E-02 


1.53E-02 


1.56E-02 


1.45E-02 


1.26E-02 


1.66e-02 


1.41e-02 


1.28e-02 


58Fe 


7.14E-06 


1.41E-05 


3.21E-05 


3.63E-05 


7.37E-05 


5.33E-05 


5.09E-03 


8.02E-05 


1.27E-07 


6.05E-05 


2.83E-05 


2.13E-04 


7.02E-06 


8.70E-06 


7.85e-05 


7.72e-05 


7.70e-05 


59C0 


4.18E-04 


4.97E-04 


5.62E-04 


5.98E-04 


6.12E-04 


5.59E-04 


9.96E-04 


5.34E-04 


1.48E-04 


5.23E-04 


5.28E-04 


5.03E-04 


6.23E-04 


6.20E-04 


5.09e-04 


4.40e-04 


7.53e-04 



58Ni 


7.26E-02 


6.78E-02 


6.98E-02 


7.13E-02 


6.63E-02 


6.89E-02 


7.54E-02 


6.90E-02 


3.81E-02 


7.01E-02 


7.29E-02 


6.26E-02 


7.48E-02 


6.16E-02 


5.90e-02 


5.09e-02 


5.01e-02 


60Ni 


2.01E-03 


2.63E-03 


3.22E-03 


3.58E-03 


4.31E-03 


4.40E-03 


8.79E-03 


4.54E-03 


7.64E-04 


4.56E-03 


4.67E-03 


4.57E-03 


5.79E-03 


6.75E-03 


4.87e-03 


5.39e-03 


5.36e-03 


6lNi 


8.45E-05 


9.36E-05 


1.02E-04 


l.OlE-04 


9.18E-05 


6.57E-05 


1.20E-04 


4.76E-05 


3.09E-05 


3.51E-05 


1.69E-05 


2.56E-05 


1.48E-05 


1.61E-05 


5.04e-05 


4.62e-05 


2.71e-05 


62Ni 


6.55E-04 


7.31E-04 


8.10E-04 


8.06E-04 


7.49E-04 


5.27E-04 


2.85E-03 


4.00E-04 


2.46E-04 


2.94E-04 


1.31E-04 


3.01E-04 


9.83E-05 


1.05E-04 


3.01e-04 


1.54e-04 


1.25e-04 


64Ni 


4.24E-09 


8.11E-09 


1.94E-08 


2.15E-08 


4.72E-08 


2.40E-08 


9.61E-06 


3.83E-08 


5.26E-14 


2.44E-08 


5.97E-09 


1.63E-07 


1.78E-11 


2.55E-11 


3.72e-08 


3.65e-08 


3.63e-08 



63Cu 


4.44E-07 


4.94E-07 


5.59E-07 


5.71E-07 


5.62E-07 


4.43E-07 


2.32E-06 


3.78E-07 


1.72E-07 


3.20E-07 


2.10E-07 


3.23E-07 


1.91E-07 


2.29E-07 


3.17e-07 


3.27e-06 


1.29e-05 


65cu 


1.09E-07 


1.15E-07 


1.32E-07 


1.32E-07 


1.31E-07 


1.05E-07 


3.23E-07 


8.02E-08 


6.04E-08 


6.14E-08 


2.58E-08 


4.13E-08 


2.41E-08 


2.77E-08 


9.98e-08 


1.26e-07 


1.45e-07 


64zn 


1.72E-06 


1.83E-06 


2.09E-06 


2.11E-06 


2.09E-06 


1.71E-06 


1.94E-06 


1.32E-06 


9.19E-07 


1.06E-06 


5.45E-07 


6.42E-07 


5.66E-07 


6.38E-07 


2.22e-06 


1.73e-05 


2.64e-05 


66Zn 


2.55E-06 


2.70E-06 


3.11E-06 


3.13E-06 


3.04E-06 


2.41E-06 


2.74E-06 


1.79E-06 


1.44E-06 


1.39E-06 


5.97E-07 


7.64E-07 


5.74E-07 


6.45E-07 


1.43e-06 


6.29e-07 


4.37e-07 


67zn 


1.48E-09 


1.53E-09 


1.77E-09 


1.79E-09 


1.75E-09 


1.42E-09 


2.26E-09 


1.09E-09 


8.79E-10 


8.61E-10 


3.91E-10 


4.91E-10 


3.76E-10 


4.24E-10 


7.08e-10 


2.47e-10 


1.73e-10 



"^Zn 9.01E-10 9.38E-10 l.lOE-09 1.12E-09 1.13E-09 8.87E-10 1.07E-08 7.24E-10 5.26E-10 5.68E-10 2.58E-10 5.91E-10 2.29E-10 2.53E-10 3.70e-10 1.71e-10 1.51e-10 

'OZn 3.39E-16 6.43E-16 1.52E-15 1.63E-15 3.67E-15 1.36E-15 2.23E-12 2.23E-15 5.97E-24 1.25E-15 1.76E-16 1.48E-14 7.40E-22 1.25E-21 2.16e-15 2.10e-15 2.09e-15 

'^Ga 5.66E-14 8.72E-14 1.77E-13 2.00E-13 3.30E-13 2.72E-13 1.04E-11 3.93E-13 6.71E-15 3.27E-13 1.39E-13 8.44E-13 9.10E-15 1.15E-14 3.81e-13 3.73e-13 3.72e-13 

'iGa 2.06E-15 2.25E-15 3.31E-15 3.49E-15 4.73E-15 3.43E-15 3.47E-13 4.03E-15 1.22E-15 3.15E-15 1.24E-15 9.47E-15 6.14E-16 7.02E-16 3.06e-15 2.69e-15 2.66e-15 
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Table 3. Nucleosynthetic yields (in solar masses) of select radioactive nuclides at time t ^ 100 s. 





Nl 


N3 


N5 


NIO 


N20 


N40 


NIOOH 


NlOO 


NIOOL 


N150 


N200 


N300C 


N1600 


N1600C 


N100_Z0.5 


N100_Z0.1 


N100_Z0.01 




2.25E-06 


7.64E-06 


6.50E-06 


2.96E-06 


6.39E-06 


3.30E-06 


3.26E-06 


2.47E-06 


3.06E-06 


1.41E-05 


1 .02E-05 


7.01E-06 


7.96E-06 


1.44E-05 


1 .20e-06 


1 .50e-07 


6.57e-10 


22 Na 


1.90E-09 


4.64E-09 


4.99E-09 


3.03E-09 


6.99E-09 


5.11E-09 


4.67E-09 


4.27E-09 


5.25E-09 


1.17E-08 


1.46E-08 


8.67E-09 


1.18E-08 


2.20E-08 


4.43e-09 


4.60e-09 


7.84e-09 


26 Al 


2.35E-07 


5.83E-07 


6.13E-07 


3.86E-07 


8.62E-07 


6.57E-07 


5.72E-07 


5.68E-07 


7.16E-07 


1.42E-06 


1.82E-06 


1.08E-06 


1.51E-06 


2.70E-06 


3.73e-07 


2.70e-07 


9.30e-08 


32Si 


6.49E-09 


1.61E-08 


1.52E-08 


8.25E-09 


1.78E-08 


1.14E-08 


9.99E-09 


9.47E-09 


l.lOE-08 


3.48E-08 


3.04E-08 


2.01E-08 


2.42E-08 


4.38E-08 


1.58e-09 


1.53e-ll 


4.58e-15 


32p 


1.46E-07 


2.55E-07 


2.99E-07 


2.65E-07 


4.80E-07 


5.01E-07 


3.77E-07 


4.96E-07 


5.88E-07 


6.68E-07 


1.03E-06 


6.29E-07 


9.56E-07 


1.42E-06 


1.89e-07 


2.43e-08 


2.02e-10 


33p 


l.OOE-07 


1.67E-07 


2.01E-07 


1.93E-07 


3.38E-07 


3.74E-07 


2.74E-07 


3.76E-07 


4.54E-07 


4.44E-07 


7.57E-07 


4.56E-07 


7.17E-07 


1.05E-06 


1.39e-07 


1.33e-08 


7.44e-ll 




1.78E-07 


2.87E-07 


3.59E-07 


3.08E-07 


5.52E-07 


5.60E-07 


4.16E-07 


5.39E-07 


6.74E-07 


7.86E-07 


1.19E-06 


7.13E-07 


1.09E-06 


1.71E-06 


8.66e-08 


3.58e-09 


2.53e-ll 


36C1 


2.15E-07 


3.24E-07 


4.13E-07 


3.98E-07 


6.72E-07 


7.66E-07 


5.61E-07 


7.77E-07 


9.58E-07 


9.13E-07 


1.48E-06 


9.15E-07 


1.43E-06 


2.07E-06 


2.81e-07 


3.59e-08 


8.47e-10 


37 Ar 


7.11E-06 


9.32E-06 


1.30E-05 


1.49E-05 


2.19E-05 


3.06E-05 


2.38E-05 


3.43E-05 


4.39E-05 


3.06E-05 


4.52E-05 


3.54E-05 


5.08E-05 


5.55E-05 


2.25e-05 


9.11e-06 


3.51e-06 


39 Ar 


4.13E-09 


5.53E-09 


7.30E-09 


6.96E-09 


1.13E-08 


1.28E-08 


9.56E-09 


1.29E-08 


1.57E-08 


1.58E-08 


2.39E-08 


1.50E-08 


2.28E-08 


3.30E-08 


2.46e-09 


1.19e-10 


1.47e-12 


40k 


1.76E-08 


2.39E-08 


3.18E-08 


3.07E-08 


4.89E-08 


5.68E-08 


4.21E-08 


5.81E-08 


7.39E-08 


6.86E-08 


1.06E-07 


6.72E-08 


1.04E-07 


1.48E-07 


1.48e-08 


1.30e-09 


5.80e-ll 


4lCa 


1.28E-06 


1.65E-06 


2.30E-06 


2.65E-06 


3.80E-06 


5.37E-06 


4.22E-06 


6.07E-06 


7.86E-06 


5.35E-06 


7.79E-06 


6.20E-06 


8.92E-06 


9.63E-06 


3.85e-06 


1.42e-06 


4.92e-07 


44Ti 


3.24E-06 


4.38E-06 


5.22E-06 


6.47E-06 


8.11E-06 


9.10E-06 


8.88E-06 


9.98E-06 


l.llE-05 


1.03E-05 


6.37E-06 


1.19E-05 


6.00E-06 


3.96E-06 


1.08e-05 


1.18e-05 


1.27e-05 


48 Y 


1.88E-08 


2.47E-08 


3.53E-08 


4.29E-08 


6.23E-08 


8.03E-08 


7.01E-08 


9.12E-08 


1.09E-07 


8.69E-08 


9.80E-08 


9.77E-08 


1.05E-07 


l.lOE-07 


5.72e-08 


3.84e-08 


3.01e-08 


49 Y 


6.80E-08 


9.28E-08 


1.31E-07 


1.63E-07 


2.41E-07 


3.15E-07 


3.21E-07 


3.57E-07 


3.93E-07 


3.51E-07 


3.72E-07 


4.00E-07 


4.00E-07 


4.01E-07 


1.95e-07 


9.93e-08 


6.32e-08 



48Cr 


1.02E-04 


1.41E-04 


1.66E-04 


2.07E-04 


2.67E-04 


2.85E-04 


2.86E-04 


3.14E-04 


3.36E-04 


3.30E-04 


1.79E-04 


3.77E-04 


1.49E-04 


9.24E-05 


3.43e-04 


3.73e-04 


3.87e-04 


49Cr 


7.99E-06 


1.06E-05 


1.27E-05 


1.59E-05 


2.04E-05 


2.32E-05 


2.29E-05 


2.56E-05 


2.72E-05 


2.66E-05 


1.68E-05 


3.01E-05 


1.58E-05 


9.89E-06 


2.05e-05 


1.26e-05 


6.76e-06 


51Cr 


1.39E-06 


2.12E-06 


3.13E-06 


3.89E-06 


5.95E-06 


8.11E-06 


1.13E-05 


9.29E-06 


7.34E-06 


9.13E-06 


l.llE-05 


1.03E-05 


1.34E-05 


1.42E-05 


5.26e-06 


3.53e-06 


3.36e-06 


51 Mn 


3.27E-05 


3.91E-05 


4.67E-05 


5.76E-05 


7.41E-05 


8.66E-05 


8.61E-05 


9.52E-05 


9.39E-05 


9.91E-05 


7.06E-05 


1.08E-04 


6.74E-05 


4.58E-05 


7.06e-05 


3.58e-05 


2.23e-05 


52 Mn 


1.44E-06 


1.70E-06 


2.10E-06 


2.63E-06 


3.80E-06 


4.57E-06 


4.82E-06 


5.18E-06 


4.77E-06 


5.45E-06 


4.58E-06 


5.87E-06 


4.60E-06 


3.95E-06 


4.61e-06 


4.29e-06 


4.14e-06 



53 Mn 


3.06E-05 


5.07E-05 


7.30E-05 


9.66E-05 


1.53E-04 


2.01E-04 


3.78E-04 


2.35E-04 


8.55E-05 


2.48E-04 


2.78E-04 


2.67E-04 


3.56E-04 


3.95E-04 


1.96e-04 


1.77e-04 


1.72e-04 


54Mn 


2.51E-07 


5.52E-07 


9.36E-07 


1.22E-06 


2.02E-06 


2.44E-06 


1.58E-05 


3.03E-06 


2.60E-07 


3.07E-06 


3.04E-06 


3.86E-06 


3.15E-06 


4.02E-06 


2.84e-06 


2.77e-06 


2.77e-06 


52Fe 


3.49E-03 


3.94E-03 


4.36E-03 


5.20E-03 


6.80E-03 


7.18E-03 


7.29E-03 


7.93E-03 


7.91E-03 


8.27E-03 


4.98E-03 


9.01E-03 


4.06E-03 


3.03E-03 


8.40e-03 


8.90e-03 


9.10e-03 


53Fe 


4.94E-04 


4.76E-04 


5.17E-04 


6.11E-04 


7.92E-04 


9.12E-04 


9.08E-04 


l.OOE-03 


9.12E-04 


1.04E-03 


7.76E-04 


1.08E-03 


6.99E-04 


5.25E-04 


8.17e-04 


5.81e-04 


4.51e-04 


55Fe 


1.33E-04 


3.40E-04 


5.33E-04 


7.40E-04 


1.20E-03 


1.59E-03 


2.94E-03 


1.86E-03 


2.75E-04 


2.01E-03 


2.37E-03 


2.12E-03 


3.08E-03 


3.29E-03 


1.73e-03 


1.66e-03 


1.66e-03 



59Fe 


2.93E-10 


5.80E-10 


1.34E-09 


1.47E-09 


3.16E-09 


1.71E-09 


5.10E-07 


2.72E-09 


1.04E-15 


1.81E-09 


5.03E-10 


1.07E-08 


5.75E-13 


9.55E-13 


2.65e-09 


2.60e-09 


2.59e-09 


60Fe 


4.95E-11 


1.03E-10 


2.45E-10 


2.67E-10 


6.17E-10 


2.60E-10 


2.43E-07 


4.20E-10 


3.64E-18 


2.47E-10 


4.91E-11 


2.35E-09 


3.02E-15 


5.33E-15 


4.07e-10 


3.98e-10 


3.96e-10 


55Co 


8.61E-03 


6.58E-03 


6.26E-03 


6.86E-03 


8.48E-03 


1.05E-02 


1.03E-02 


1.14E-02 


9.01E-03 


1.19E-02 


1.08E-02 


1.13E-02 


l.OOE-02 


8.03E-03 


9.33e-03 


7.04e-03 


6.17e-03 


56Co 


5.08E-05 


4.84E-05 


5.45E-05 


6.44E-05 


8.25E-05 


1.08E-04 


1.22E-04 


1.18E-04 


5.43E-05 


1.26E-04 


1.39E-04 


1.19E-04 


1.48E-04 


1.21E-04 


1.07e-04 


9.92e-05 


9.66e-05 


57Co 


5.71E-05 


1.58E-04 


2.50E-04 


3.52E-04 


5.69E-04 


7.44E-04 


1.36E-03 


8.70E-04 


6.54E-05 


9.52E-04 


l.llE-03 


9.92E-04 


1.42E-03 


1.50E-03 


8.43e-04 


8.28e-04 


8.25e-04 



58Co 


3.11E-07 


7.94E-07 


1.30E-06 


1.79E-06 


2.88E-06 


3.63E-06 


1.07E-05 


4.35E-06 


1.23E-07 


4.72E-06 


5.04E-06 


5.14E-06 


6.08E-06 


7.29E-06 


4.27e-06 


4.22e-06 


4.21e-06 


60Co 


2.01E-09 


3.77E-09 


8.36E-09 


9.40E-09 


1.78E-08 


1.36E-08 


8.78E-07 


2.03E-08 


6.97E-13 


1.59E-08 


6.87E-09 


4.97E-08 


2.04E-10 


3.00E-10 


1.99e-08 


1.96e-08 


1.96e-08 


56Ni 


l.llE-00 


1.04E-00 


9.74E-01 


9.39E-01 


7.78E-01 


6.55E-01 


6.94E-01 


6.04E-01 


5.32E-01 


5.66E-01 


4.15E-01 


5.12E-01 


3.64E-01 


3.22E-01 


6.29e-01 


6.49e-01 


6.55e-01 


57Ni 


3.36E-02 


3.04E-02 


2.89E-02 


2.78E-02 


2.25E-02 


1.96E-02 


2.06E-02 


1.79E-02 


1.42E-02 


1.69E-02 


1.42E-02 


1.45E-02 


1.31E-02 


l.llE-02 


1.57e-02 


1.33e-02 


1.19e-02 


59Ni 


8.22E-05 


1.32E-04 


1.78E-04 


2.22E-04 


3.01E-04 


3.57E-04 


5.67E-04 


3.93E-04 


5.78E-05 


4.22E-04 


4.80E-04 


4.23E-04 


5.82E-04 


5.77E-04 


3.80e-04 


3.63e-04 


4.14e-04 


63Ni 


1.85E-09 


3.43E-09 


7.79E-09 


8.72E-09 


1.74E-08 


1.13E-08 


1.38E-06 


1.76E-08 


1.58E-14 


1.27E-08 


4.26E-09 


5.31E-08 


1.85E-11 


2.92E-11 


1.72e-08 


1.69e-08 


1.68e-08 


62zn 


6.48E-04 


7.18E-04 


7.80E-04 


7.72E-04 


6.86E-04 


4.74E-04 


5.25E-04 


3.22E-04 


2.45E-04 


2.31E-04 


9.92E-05 


1.34E-04 


9.39E-05 


9.86E-05 


2.24e-04 


7.89e-05 


5.01e-05 


65zn 


3.56E-10 


4.88E-10 


6.69E-10 


7.18E-10 


8.76E-10 


8.30E-10 


1.62E-09 


7.35E-10 


4.09E-10 


6.35E-10 


4.11E-10 


5.23E-10 


3.95E-10 


5.08E-10 


8.46e-10 


9.95e-10 


l.lle-09 


65Ge 


5.98E-08 


6.80E-08 


8.19E-08 


8.30E-08 


8.57E-08 


7.07E-08 


7.22E-08 


5.29E-08 


4.32E-08 


3.99E-08 


1.66E-08 


2.07E-08 


1.69E-08 


2.03E-08 


6.64e-08 


8.45e-08 


9.79e-08 


68Ge 


8.91E-10 


9.19E-10 


1.06E-09 


1.07E-09 


1.03E-09 


8.30E-10 


8.90E-10 


6.33E-10 


5.26E-10 


5.05E-10 


2.40E-10 


2.84E-10 


2.29E-10 


2.53E-10 


2.81e-10 


8.32e-ll 


6.38e-ll 



EXC 153. FKR, MF and SAS acknowledge travel support by the 
DAAD/G08 German- Australian exchange program. 
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